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We present an analytical method which enables one to find the exact spatial dependence of the 
indirect RKKY interaction between the localized moments via the conduction electrons for the 
arbitrary dimensionality n. The corresponding momentum dependence of the Lindhard function is 
exactly found for any n as well. Demonstrating the capability of the method we find the RKKY 
interaction in a system of metallic layers weakly hybridized to each other. Along with usual 2/cf in- 
plane oscillations the RKKY interaction has the sign-reversal character in a direction perpendicular 
to layers, thus favoring the antiferromagnetic type of layers' stacking. 
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The Ruderman-Kittel-Kasuya-Yosida (RKKY) inter- 
action was found to play an important role in various 
problems involving the interaction of the localized mo- 
ments in a metal via polarization of conduction electrons. 
In this paper we provide the exact derivation of the spa- 
tial and momentum dependence of RKKY interaction for 
arbitrary dimensionality. In three dimensions the sp 
tial dependence of this interaction was obtained in Ref. 
about forty years ago. It took some time to obtain the 
corresponding exact expression in two dimensions .□ The 
primary goal of this paper is thus to present a promis- 
ing analytical method for evaluation of the corresponding 
expressions. 

The obtained expressions being the analytical function 
of dimensionality might prove to be useful for theoreti- 
cal approaches dealing with the fractional and infinitely 
large dimensions. Closer to practical needs, we apply 
our method for an analysis of multi-layer metal struc- 
ture. We demonstrate here the sign-reversal character 
of the interaction in a direction perpendicular to layers. 
The period of these oscillations coincides with a double 
interlayer spacing thus favoring the antiferromagnetic or- 
dering of layers. It is interesting to note that the above 
alternation of sign of the interaction coexists with the 
usual in-plane 2k p oscillations. 



I. RKKY INTERACTION IN N DIMENSIONS 

We begin with conventional form of the exchange in- 
teraction between the localized moment J and electron 
spin density s(r) : 

V(r) = -AJ(R)s(r)5(R-r) (1) 

Here A is the exchange coupling constant. The RKKY 
interaction between two localized moments via the con- 
ducting electrons may then be written in the following 
form 



H, 



HKKY — — 2^ 2 JlJ2X(R'12) 



(2) 



where R— dependent part of the interaction coincides 
with the Fourier transform of the non-uniform static sus- 
ceptibility x{l) (Lindhard function) and is given by : 



X(R) = -T £ <?(«,;, ,R) 



(3) 



here Matsubara frequency u>i = nT(2l + 1) and the elec- 
tronic Green function is 



G(iw,R) = 



d n k exp(ikR) 

(2n) n iuo — £fe 



(4) 



We focus our attention below at the case of low 
temperatures and use the limiting relation — > 

The quadratic electron dispersion in n dimensions is 
explicitly assumed : 



£k = k 2 /2m — /i 



(5) 



with the Fermi energy \i = k F /2m. First we use the 
following representation of the Green function : 



oo 



ikR 



2m 



(6) 



where we introduced the value z = fj, + iu> and a — 
sgn(uj)'K/2. The Gaussian integration over k gives 



G(iu,R) = (£) n/2 e -«( 1 +™/ 2 ) 



dr 

7-»/2 



exp|rze 



(7) 



with p — mR 2 . We notice that the last integral can bs 
expressed via the modified Bessel (Mcdonald) function,E 
namely 



m \ "+ 1 / \J~2zp 



) [^J^) K v (yp2zp), (8) 



we defined v = n/2 — 1 here. In this equation the 
branch of root \J—2zp should be chosen from the con- 
dition of its positive real part. In particular, this latter 
condition means that the argument of Mcdonald function 
K V (\J— 2zp) has a discontinuity at u = : 



1 



K v {^2z~p) = 



e- um / 2 Hi 2) {k F R) 



to 



+0 
-0 



where H^' 2 \x) are Hankel functional. 

Next we observe that one can change the variable to — > 
Z = y / —2zp in (^) and integrate over complex Z using 
the exact form for the Green function (H). Note that 
the limits of integration by Z are (1 ± i)oo. Without 
the above discontinuity at Z = ±iy/2jlp one could shift 
the integration contour to Z — » +oo and obtain zero for 
(||) in view of the property K V {Z) oc e~ z . Due to the 
discontinuity, the function \ has a finite value. After 
some calculations we gcta 



x(R) 



( _kF_Y l d2 
2nR J 



R 2 $ n (k F R) 



(9a) 



$ n {x) = J n /2-i(x)Y n/2 _ 1 (x) + J n/2 (x)Y n/2 (x) (9b) 

This expression is the main finding of this Section. Let 
us take a closer look at this result. First we note that 
(^) is the continuous function of both distance R and the 
dimensionality n. 

At large distances kpR ^> 1 the leading terms of the 
asymptotes of Bessel functions appearing in (j^) cancel 
each other. Next terms produce the following expression^ 



x(R) 



i(2^)" sin(2fc ^ + ™ /2) (l0) 



In particular cases of physical interest the general 
expression (^) immediately provides the exact iorm of 
RKKY interaction in three and two dimensionsEra. For 
n = 3 one has 



x(R) 



mkp ( , sin2fcpi? N . 

' cos2k F R * ) , (11) 



87T 3 i? 3 



2k F R 



and for n = 2 

X (R) = -^[Jo(k F R)Y {k F R) + J 1 (k F R)Y 1 (k F R)]. 

47T 

(12) 

The one-dimensional case can be obtained either by the 
continuationo of Bessel functions in (^) upon the index 
n or by the direct evaluation of the, integral with Green 
function (ph at n = 1. The result isB 



X (R) = -si(2k F R), (n = l) 
with the sine integral 

si(a;) = 



(13) 



00 dt 
— sint 



t 



It is useful to define here the density of states 
at the Fermi level N(E F ) = J d n k/ (2Tr) n S(e k ) = 
—ir~ 1 ImG(iuj — > +iO,R — ► 0). From (^) one immedi- 
ately finds 



N(E F ) 



27rr[r7,/2] \Att 



(14) 



Now knowing the exact expression (g) for RKKY inter- 
action in R— space one can find its correspondence in 
q— space as follows. 

X (q) = J d n Ke^ R X (R) (15) 

We see that x(l) is reduced to the Mellin convolution 
of Jv{x) and $„(x). A straightforward calculation gives 
then the answer expressed via the Gauss hypergeometric 
function 2-Fi[a, b, c; z\. 

(16a) 



X (q) = N(E F )cf> n ( ^~ 



F 



(t>n{x) 



1 

X '2 ;1 



n 

r 



F 
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x > 1 

X < 1 



(16b) 



Again, the result is the continuous function both in q 
and n. From the general properties of hypergeometric 
function, one has ^>„(0) = 1, (f> n (x » 1) ~ l/(nx^ and 
</>„(l) = l/(n — 1). At last, one easily verifyd that 
the expressions known previouslyH'Bo are reproduced in 
particular cases n = 1, 2, 3. 



II. RKKY INTERACTION FOR QUASI-2D CASE 

It is interesting to note that both exact (||) and asymp- 
totic ([h]) expressions for x(R) let one mimic the "switch- 
ing on" the extra dimensionality of a metal by simple 
change of the index n. Thus at first sight one could tackle 
the case of system of weakly hybridized metallic planes 
by ascribing the dimensionality 2 + e to it. Actually the 
situation is more complicated as we discuss below. 

Let us consider the (infinite) set of metallic layers, 
weakly connected to each other. By this we assume the 
following dispersion : 



£k = {k 2 x + k 2 y )/2m - (i - C cos k z 



(17) 



with £ -C n and \k z \ < ir . The Fermi surface has 
a cylinder-like shape with maximum and minimum in- 
plane radii defined by kp = y^2m(fi ± £). We write 
kR = k||R|| + k z l where I is the integer number of lay- 
ers. Below we retain the definitions of p = 77ii?| and 
k F = y / 2ro7t for the simplicity of writing. 



Using (||), (0) we come to expression 
G(R) 



in 



-ial 



2tt 

00 dr 



Ji(T()exp[Tze l 



o 



2t 



(18) 
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The latter integral can be evaluated ( for large in-plane 
distances ) by the steepest descent method. We note 
that when \J2zp ~ fepiZii 3> 1 , the principal contri- 
bution to the integral comes from the vicinity of the 
point To = y/ p/2z ~ kpRn/2/j,, more rigorously, at 
t = r (l + (k F i?||)~ 1/2 0(l)). It follows then that at 
k F R\\ ^ (/VC) 2 one can replace r by r in an argument 
of the Bessel function in (jig ) . As a result the quasi- two- 
dimensional RKKY interaction is 

X(R) = X2d(R\\) J? (R\\/Ro) (-1)', R\\ < k F R 2 (19) 



where k F Ro = 2/i/C ^ 1 and X2d(R) is given by (|12[). 
Analyzing this expression we first note the appearance of 
the length scale Rq inversely proportional to the strength 
of hybridization of layers. Since Ji(x) ~ (x/2) /l\ for 
x < 1, at moderate in-plane distances Ru < i?o the in- 
teraction ( |l9| ) rapidly decays as a function of I, 

More interesting however is the fact of sign-reversal 
character of interaction in a direction perpendicular 
to layers. This modulation of interaction has a period 
exactly coinciding with a double lattice parameter and 
should obviously lead to the preferential antiferromag- 
netic stacking of layers. This phenomenon is accom- 
panied by the usual 2k F oscillations of in-plane term 

X2d{R\\)- 

We can further clarify this point by performing the 
Fourier transform with the resultlx 



X (k) = / d 2 R||e ik ll R ii X2D (i*||) Jo ( 2^cos^ 



R 



Rq 



(20) 



with fc* = (2/Rq) cos k z /2 <C 1. The last integration is 
over the angle ip of k* in the plane, i.e. the points k|| +k» 
lie on the circle of radius fc* and with center at k|| . 

According to (j^) a maximum of X2Z>(kii) at some ko, 
corresponds to the possible in-plane magnetic ordering, 
characterized by this wave-vector. We see from ( ^o|) that 
the inclusion of weak interplane hopping leads to the po- 
sition of true maximum at k z = 7r, when fc* = 0. The 
other values of k z cause the loss in the magnetic energy 
of order of the value x(ko)fc 2 /fcf- 

It is worth noting that the weak interplane hopping 
£ and large effective in-plane Fermi momentum k F are 
obviously realized in the high-T c cuprates. It is known 
that the 2D. Fermi surface has a complicated form in these 
substancesB. We believe however that RKKY interaction 
in this case preserves the general form (|l9|), j20| ) with 
anisotropic in-plane form of interaction X2u(kii) and Rq 
defined by some effective k F ~ 1. It is also known that 
in the compounds RBa2Cu307_5 the subsystem of rare- 
earth ions undergoes a magnetic ordering transition at 
low temperaturesQ ; generally the type of ordering de- 
pends on a particular ion R 3+ . Remarkable fact is how- 
ever that for all substances the antiferromagnetic stack- 
ing of magnetic R 3+ layers was reported in accordance 
with our finding (n9J) . 



Concluding this Section, we wish to stress the following 
point. It was previously showrolij for the case of compli- 
cated Fermi surface (FS), that the period of oscillations 
( and general power-law behavior ) of the RKKY inter- 
action is determined by the calipering pairs of points on 
the FS. These are the points where the direction of nor- 
mal to the FS is (anti)parallel to the direction of R. One 
can see that the very notion of calipering points implies 
the closeness of the Fermi surface at a given direction of 
R. On the contrast, the FS is obviously open in our case 
at the z direction and the oscillations exist, albeit the 
roughly exponential law of their decay. 

In conclusion, we found the exact form of spatial de- 
pendence of RKKY interaction for arbitrary dimension- 
ality. Its counterpart in momentum space is also found. 
Applying our method to the system of weakly hybridized 
metallic layers we demonstrate the existence of spatial 
oscillations of indirect RKKY exchange in the direction 
perpendicular to layers. The period of oscillations equals 
exactly to double interlayer spacing, which indicates the 
preferential antiferromagnetic ordering of layers. 
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